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Abstract 

We compute the 1-loop (a 2 ) correction to hard spectator scattering in non-leptonic 
B decay tree amplitudes. This forms part of the NNLO contribution to the QCD 
factorization formula for hadronic B decays, and introduces a new rescattering 
phase that corrects the leading-order result for direct CP asymmetries. Among 
the technical issues, we discuss the cancellation of infrared divergences, and the 
treatment of evanescent four-quark operators. The infrared finiteness of our result 
establishes factorization of spectator scattering at the 1-loop order. Depending on 
the values of hadronic input parameters, the new 1-loop correction may have a 
significant impact on tree-dominated decays such as B —> tttt. 


1 Introduction 


The majority of observables at the B factories is connected with branching fractions 
and CP asymmetries of hadronic B decays to two charmless mesons, for which strong- 
interaction effects are essential. There is some control over these effects, since the decay 
amplitudes factorize in the heavy-quark limit. In the QCD factorization framework [L 
the matrix elements of the effective weak interaction operators take the (schematic) 
expression 

(M 1 M 2 \Qi\B) = F BMl (0) T] * /m 2 0m 2 + Tf * Mb+ * * /m 2 0m 2 - (1) 

The long-distance strong-interaction effects are now confined to a form factor F BMl (0) 
at q 2 = 0, decay constants /m, and light-cone distribution amplitudes 0 m- The benefit 
is that information extraneous to two-body B decays is available for these, and that 
the short-distance kernels Zy 11 can be expanded in a perturbation series in the strong 
coupling a s . Both kernels are currently known from Tj at order a s . While for T- 
this includes a 1-loop correction to “naive factorization”, in case of T} 1 the order a s 
contribution is actually the leading term. It originates from the tree-level exchange of a 
hard-collinear gluon with the spectator-quark in the B meson, as indicated in Figure |2] 
below. (The class of corrections from fermion-loop insertions into the gluon propagator 
is also known [2j. In spectator scattering these n/-terms are all connected with the 
hard-collinear scale and make no contribution to the hard 1-loop correction, which we 
compute here.) 

In this paper we shall compute the 1-loop (a^) correction to the spectator-scattering 
kernel T- ] for what is known as the (topological) “tree amplitudes” in two-body decays. 
There are several motivations for performing this calculation: 

• As in any perturbative QCD calculation the 1-loop correction is necessary to elim¬ 
inate scale ambiguities. In the present case of spectator scattering the characteris¬ 
tic scales are mi, and (mi, Aq CD ) 1//2 . The latter being only about 1.5 GeV, a 1-loop 
calculation is necessary to ascertain the validity of a perturbative treatment by 
showing that the expansion converges. The 1-loop correction to spectator scat¬ 
tering forms part of the next-to-next-to-leading order (NNLO) contribution to the 
decay amplitudes. 

• At order a s the strong interaction phases, and hence direct CP asymmetries, orig¬ 
inate entirely from the imaginary part of the kernel Tj in the first term on the 
right-hand side of Q. The 1-loop correction to T} 1 introduces a new rescatter¬ 
ing mechanism by spectator scattering. Its calculation represents an important, 
presumably dominant, part of the next-to-leading order (NLO) result for the CP 
asymmetries. The NLO result will be needed to resolve or understand potential 
discrepancies of the LO result with experimental data. 

The organization of the paper is as follows. In Section [21 we set up the definitions and 
matching equations for the calculation of the hard-scattering kernel Tj 1 , which is then 
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Figure 1: Flavour and colour flow for insertions of Q\. The connected fermion lines 
indicate the contraction of spinor indices. Left: “right insertion”. Right: “wrong inser¬ 
tion”. 

described in Section 01 The expression for the kernel is given at the end of that section. 
In Section 01 we obtain the tree amplitudes a\ i2 {MiM 2 ) in a convenient representation, 
where the light-cone distribution amplitudes are integrated in the Gegenbauer expansion. 
The numerical effect of the new correction on the tree amplitudes and the B —> hit 
branching fractions is investigated in Section 01 We conclude in Section 01 

2 Set-up and matching 

2.1 Flavour and colour 

We are concerned with the current-current operators in the effective weak Hamiltonian 
for b —> u transitions given by 

'HeS — VuD^ub {C\Q\ + C 2 Q 2 ) + h.c., 

Ql = [UaYil - 75- 75 )u b \, 

Q 2 = [UbYi 1 - lb)b a }[Dal^{l - 75 M, (2) 

with a, b denoting color, and D = d or s. There are two possible flavour flows to the final 
state as illustrated in Figure[T]for Q\. In case of the colour-allowed tree amplitude (left), 
denoted ai(MiM 2 ) following the notation 1 of 21, meson M 2 represented by the up-going 
quark lines has the flavour quantum numbers of [uD], and Mi those of [q s u], where q s 
denotes the flavour of the spectator anti-quark in the B meson. In case of the colour- 
suppressed tree amplitude a 2 {MiM 2 ) (right in Figure [TJ), d ie corresponding quantum 
numbers are [uu] and [q s D] , respectively. In addition there exist “penguin contractions”, 
where the u and u fields from Qi are contracted in the same fermion loop. Together 
with other operators from the effective Hamiltonian they contribute to the (topological) 
penguin amplitudes, which we do not consider in this paper. Thus, in the computation 

1 The normalization is such that at tree-level a\{MiM 2 ) = C\ + C 2/3 and U 2 (MiM 2 ) = C 2 + Ci/3 
with Ci ~ 1.1 and C 2 ~ —0.2. See section 2.2, for the relation between decay amplitudes and the 
cu parameters. 
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of the (topological) tree amplitudes there appear four short-distance coefficients, two 
corresponding to the matrix element of Qi as shown in Figure 0 and two corresponding 
to Q 2 , which differ only by the colour labels at the operator vertex. It will be seen from 
the final result that only two of the four coefficients are different, because Q\ and Q 2 are 
equivalent by a Fierz transformation, when the flavours u and D are not distinguished. 
However, since we use dimensional regularization, Fierz symmetry cannot be assumed 
to hold a priori 

We shall refer to the flavour-flow diagrams, where the spinor indices are contracted 
along the quark lines of M 2 (left diagram of Figure^]), as the “right insertions” of Q 12 ; the 
other contraction (right diagram) is the “wrong insertion”. Exactly the same diagrams 
contribute to the two right (wrong) insertions, only the colour factor is different for 
each diagram, since the two operators Q lj2 have different colour-orderings. With colour 
and flavour thus understood, we will omit colour and flavour labels in the subsequent 
discussion of operator matching. 

2.2 Matching onto SCETi 

The short-distance kernels T) I,n can be determined by extracting the hard and hard- 
collinear momentum regions from quark decay amplitudes according to the strategy of 
expanding Feynman diagrams by regions |5J. The calculation becomes more transparent, 
when it is organized as an operator matching calculation in soft-collinear effective theory 
(SCET) 0. The spectator-scattering kernel is obtained by the matching sequence QCD 
—> SCETi —> SCETn, by which hard fluctuations (k ~ mb, virtuality m 2 ) and hard- 
collinear fluctuations ( n + k ~ m b , k± ~ (m^AgcD) 1 / 2 , ~ Aqcd, virtuality m b Aqcd) 

are integrated out in two steps. This method has by now been worked out completely 
for heavy-to-light form factors at large recoil energy of the light meson, both to all 
orders urn, and by explicit 1-loop calculations of the short-distance coefficients El 13 EDI. 
For application to non-leptonic decays the effective theory has to be extended to include 
two sets of collincar fields corresponding to the (nearly) light-like directions of the two 
final-state mesons. As explained in EH this is a relatively minor complication, because 
the collincar fields for different directions decouple already at the scale mb- 

Our SCET conventions follow those of the form-factor calculations El El EDI. Meson 
Mi, which picks up the spectator anti-quark from the B meson, moves into the direction 
of the light-like vector n_. The collinear quark field for this direction is denoted by £ with 
t/_£ = 0, the corresponding collinear gluon field is A c \. The second meson moves into the 
opposite direction n + , and the collinear fields for this direction are y, satisfying rj + x = 0, 
and A c2 . The heavy quark field h v is labeled by the time-like vector v = (n_ + n + )/ 2 
with v 2 — 1. 

In [3 a power-counting argument has been developed to identify the SCETi oper¬ 
ators that can appear at leading power in the 1/mb expansion of heavy-to-light form 
factors. Applying this argument to the two collinear directions separately, we find that 
Q 12 can match to only two operators in SCETi with non-vanishing matrix elements 
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(MiM 2 | . . . | B) for non-singlet mesons Mp 2 . 2 The leading operator in the collinear-2 sec¬ 
tor is uniquely given by (xW c 2 )(£n_)j/_(l — 75 )(VEc 2 X)( 0 ). The two operators are then 
constructed by multiplying this with an AO- and a B 1-type current for the B —> Mi 
transition Due to chirality conservation and the requirement that the operator be a 
Lorentz scalar, there is only one current of each type. The two SCETi operators thus 
obtained can be arranged to reproduce the structure of the factorization formula (JTJ) by 
defining 


o\t) = (,\iv2)(in-)A(i- 75 )(iy 2 x) [cjy 1 - ls )h 


O ( t,s ) = 


1 r 

m b L 


- [ dsCf^is) m cl )i + [Wliiip Lc iW c i}(sn + ){ 1 + 7 5 K 

777 /^ J 

(xWa)(tn-)- A(1 - 7 5 )(Wlx) 


X 


m cl )^±[W^iip ±cl W cl ](sn + )( 1 + 75 )K 


( 3 ) 


The first operator includes the short-distance coefficients CCj^' 1 (s) such that its 
matrix element is proportional to the form factor f+ M ' ( 0 ) (Aq A!i ( 0 ) for vector mesons) 
in QCD (not SCETi). The expressions for the coefficients to 1-loop (more precisely, 
their momentum space Fourier transforms) can be found in [TO] . but they will not be 
needed here. In © fields without position argument are at x — 0 , and the field products 
within the large brackets are colour-singlets. We do not consider colour-octet operators, 
since their matrix elements between meson states vanish. Although the second operator 
carries an apparent l/m b suppression, both operators are in fact leading, because the 
matrix element of O z (t) is suppressed. Hence, at leading order in l/m bl the operators 
Q 1,2 from © are represented in SCETj by the equation 


Q= jdtf\i)0\t) + JdidsH n (t , s)O u (t,s) 


( 4 ) 


with s = n + p's = tubs, t = n_qt = mst, and 7 / (q) the momentum of Mi (M 2 ). Of 
the two matching coefficients T l {u) = f dte lut T l (t) is already known to the 1-loop order 
(^s) ra- In this paper we compute the 1-loop (a^) correction to 

H u (u,v) = j dtdse^+^-vW H u (t,s). (5) 

We recall that on accounting for flavour there are actually two copies of O l {t , s), O n (t, s) 
with different flavour structure, and given the two operators Qi i2 in the effective Hamil¬ 
tonian, there are four different coefficient functions H n (u, v), which we do not distinguish 
here to simplify the notation. 

2 Recall that we are not counting flavour degrees of freedom. Mesons with flavour-singlet components 
require additional two-gluon operators, as well as a term that does not factorize in SCET m- 
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To see how CD follows and to make the overall factors explicit, we evaluate the matrix 
element of CD for the case that Mi and M 2 are both pseudoscalar mesons. The SCET 
Lagrangian contains no leading-power interactions between the collinear-2 and collinear- 
1 fields after decoupling soft gluons from the collincar-2 sector by a field redefinition 
(second paper of |5|). The matrix elements of O l (t, s), O n (t,s ) fall apart into ([TDj. 
eqs. (18,81) with E = n + p '/2 = wb/2) 

(M 2 |(xir c2 )(to_)?h(l-75)(W c t 2 n|0) = f due™ 1 

A A Jo 

(^ 1 |(^^c 1 )-^[hEj' 1 *^j_ c hl / cl ](sn + )(l + ^ 5 )h v \B) = -m b m B dr e iT ° S Ml (r), (6) 
2 jo 

such that 


(M 1 M 2 \Q\B) 


im 2 B 


/+ Ml (0) r duT\u) f Ma <t> Ma (u) 

Jo 


1 

2 



dudz H ll (u , z) Smi(1 


z) /m 2 0m 2 (m) 


(7) 


A demonstration of factorization should provide an argument for the convergence of the 
various convolution integrals, an issue that is not solved to all orders in perturbation 
theory for the second term (spectator scattering) in the bracket. The convergence will 
be explicitly checked at 1-loop in our calculation. At the 1-loop order it is also easy to 
see by diagrammatic analysis that no operators other than O l {t ), O n (f, s ) are needed to 
reproduce the hard momentum regions. In particular any diagrams that match directly 
onto six-quark operators already in SCETj are power-suppressed. 


2.3 Matching onto SCETp 

To complete the derivation of (GJ the hard-collinear scale is integrated out by matching 
onto SCETn- Hard-collinear momentum regions appear only in spectator scattering, 
since an external soft momentum is required. Thus the first term in the bracket of m 
is left unchanged, while the SCETj form factor S jWl related to the matrix element Q 
must be matched onto SCETn. No new calculation is needed for this step, since we can 
use ([HU, eq. (86)) 

s Mi(r) = -— / du dvJ\\(T-,v,u)f B <l>B+(w)fM 1 <f>M 1 (v), (8) 

4 m b jo jo 

where the “jet function” J» (r; v, u>) has been calculated to 1-loop in Bum, and f B is 
^m B times the B decay constant in the static limit ([IHI, eq. (83)). Inserting this into 
0 , we obtain 

(M 1 M 2 \Q\B) = im\ | f+ Ml ( 0) duT l {u) /m 2 0m 2 (m) 
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Figure 2: Tree diagrams for H ll ^°\u,z). 


+ 


du / dudv T n (u;, u, v ) f B <t> B +(u) (v) f M2 <$>M 2 


(9) 


which (up to a normalization factor im 2 B ) is (£Q) with 


T n (u;, u, u) =- — [ dz H n (u, z) J\\(l — z;v,u). 

8m h Jo 


( 10 ) 


The jet function is unique, i.e. all four hard-scattering functions H u (u, z ) are convoluted 
with the same Jy(l — z;v,u). 

The tree-level expressions for the hard coefficient functions (when not zero) and the 
jet function are 


H im (u,z ) 



1 2 

iVc^’ 

4:7ra s CF 1 
N c m B uJV 


( 11 ) 


where we introduced the QCD colour factors C B = {N 2 — 1)/(2 N c ), Ca = N c = 3, and 
the “bar notation”, in which x = 1—x for convolution variables x. Only the two diagrams 
shown in Figure [21 have to be computed to obtain H ll ^\u, z). The other two diagrams 
with attachments to the horizontal quark lines are included in the tree contribution to 
T l Cf\ and thus belong to the T l O l term in (j3J). Combining the tree coefficients, we 


obtain 


T im (uj,u,v) 


na s CF 1 
N 2 rribujuv ’ 


( 12 ) 


which reproduces the result from jTJ. Note that m b denotes the 6-quark pole mass, and 
m B the B meson mass, but that factors of m b and m B have not been distinguished in Tj, 
since the difference is a power correction. 


3 1-loop calculation 

In this section we describe technical aspects of the computation of H ll (u , v). We calculate 
the 5-point b —> g C 2 <?c 2 Q , cifi , ci amplitude 

(q(qi)q(q 2 )q(p'i)g(p' 2 )\Q\Kp)) ( 13 ) 
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and the corresponding SCETi matrix elements of the right-hand side of (QJ). With the 
exception of one class of diagrams to be discussed below, the parton momenta can be re¬ 
stricted to their leading components. Thus for the partons in the collinear-2 direction we 
put q\ = umbn. |_/ 2 , g 2 = umpn+j 2 , for those in the collinear -1 direction p\ = vmpri-/2, 
p' 2 = hmf,n_/2, and for the heavy quark momentum p M = myiV. For such external 
momenta the SCET and HQET spinors coincide with the QCD ones. 

We use dimensional regularization with d = 4 — 2e and an anti-commuting 75 (NDR 
scheme). The amplitude (Td has ultraviolet (UV) and infrared (IR) singularities. The 
former must be subtracted in accordance with the definition of the operators Q in the 
effective Hamiltonian uni; the latter in accordance with the definition of the jet function 
J || and light-cone distribution amplitudes. This is accomplished by using MS subtractions 
and a certain prescription for dealing with evanescent operators. We first discuss the 
“right insertion” of Q, in which the quark spinor indices are contracted according to 
[xx] [£h„]. The “wrong insertion” leads to [£x] [x^v] , which differs from the desired order 
© by a Fierz transformation. 

3.1 Evanescent operators 

The calculation in dimensional regularization is complicated by the presence of evanes¬ 
cent products of Dirac matrices (products that vanish in four dimensions). When such 
products multiply 1/e poles they need special treatment. In our calculation there are 
evanescent products that multiply UV singularities. Their definition is related to the 
renormalization convention for Q. The NDR-MS scheme corresponds to setting 

'f'fYO- ~ 75) ® 7m7v7p(1 - 75 ) = (16 - 4e) 7 M (1 - 75 ) ® 7a*(! - 7s), (14) 

whenever the left-hand side multiplies an UV pole. All other products multiplying UV 
singularities can be reduced to (ED by permutations. 

The evanescent products that multiply IR poles are more complicated. Their treat¬ 
ment is related to the definition of evanescent operators of the O u (t, s) type in SCETi. 
To reduce the notation to the essentials, we strip off all the fields from O n (t, s) and 
represent it only by its Dirac structure, 

0 U (t , s)^y( l-7s) ® ^(! - 7 s)7jl> (15) 

where 

« i+ M . i- M , p na s 

7 = 2 ~ + T n + + 7 ^' 

In our 1-loop calculation we encounter the four operators 
Oi = ^(1 - 7s) ® ^(1 - 75)7 jl. 

°2 = ^7x72(1 - 75) ® ^(1 - l5h±a, 
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0 3 = y7±7 ^l(! - 75 ) <8 y (1 - 7s)7±7±a7±/3, 

0 4 = y7±7±7±7±(l - 75 ) ® y (1 - 75)7±«7 ±/j7±t- ( 17 ) 

In this notation 0\ equals O n (t, s ). One easily checks that the other three operators are 
evanescent, i.e. vanish in four dimensions. These operators will disappear from the final 
result , since we shall renormalize them such that their IR-finite matrix elements vanish, 
but they must be kept in intermediate steps, hence the matching equation 0 has to be 
extended to include all four operators on the right-hand side. 

Evanescent operators appear already at tree level. In this approximation the matrix 
element m is given by 


«2>„, = 7 (| (0l ) - -4<o 2 >) . (is) 

JS C \u uu / 

(The “right insertion” of Q 1 vanishes at tree level, because the colour-trace is zero.) The 
subscript “nf” (for “non-factorizable”) means that the “factorizable” terms that belong 
to T l O l are omitted, and only the two diagrams in Figure El are included. While one can 
simply set (O 2 ) = 0 here to recover (TTT1) . since no 1/e poles are present at tree level, the 
appearance of an evanescent operator at tree level implies that one must compute the 
mixing of O 2 into 0\ in the 1-loop calculation. 

3.2 UV renormalized 1-loop amplitude 

The calculation of the 1-loop correction to H ll (u,v) involves the diagrams shown in 
FigureEl The two lines directed upward represent the quark (anti-quark) with collinear-2 
momentum proportional to n + . The horizontal lines describe an incoming bottom quark, 
and an outgoing collinear-1 quark with momentum proportional to n_. The momentum 
of the external gluon is also in the n_ direction. The calculation of diagrams with no 
gluon lines that connect the two upper lines to the horizontal lines is not necessary, since 
the definition of O l is chosen such that these diagrams contribute only to T l {u). 

The calculation of the diagrams uses standard methods. The massive box integrals in 
dimensional regularization can be evaluated adapting the method of |14j . Alternatively, 
they can be reduced to vertex integrals, because all external momenta are linear combi¬ 
nations of only two vectors n_, n + . This observation also simplifies the tensor reduction, 
since one can use 

klkl - jRaftl. (19) 

The classes A, B of 1-particle reducible diagrams must be included in the amplitude 
calculation. The heavy-quark propagator to the right of the external gluon line in class 
A is off-shell by an amount of order m^, hence these diagrams contribute entirely to the 
short-distance coefficient. Class B is more complicated, since the light-quark propagator 
with momentum p' = p\ + p' 2 has small virtuality, hence the diagram is not completely 





Figure 3: 1-loop diagrams for (Q) n f. 
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short-distance. The non-local, long-distance contributions cancel in the matching re¬ 
lation against time-ordered products of O 1 and the SCET interaction Lagrangian as 
discussed in [5j. The local contribution to the short-distance coefficient can be extracted 
via the substitution 

if i t/i + 
p' 2 n + p' 2 

A short-cut to this conclusion is obtained by observing that we can put p[ ± = p' 2L = 
p' ± = 0 , since we do not match operators with transverse derivatives, and keep only 
= n + p' n ^/2 + n_p' n\j 2. For p' L = 0 all relevant interaction terms from the SCET 
Lagrangian vanish, hence the class B diagrams are purely short-distance. Indeed, since 
the rf_ term in the propagator does not contribute owing to the on-shell spinor to the 
right, the substitution m becomes an identity for p' L = 0 . 

Ultraviolet renormalization of the amplitude involves standard counterterms from the 
QCD Lagrangian as well as the counterterms for Q. The UV-renormalized amplitude is 
written as 



<q) ». = e (Af> + 4 1 ’) ( 21 ) 

i =1 

where (Ojf 0 -* denotes the partonic tree-level matrix element of Oj, equal to the Dirac 
matrix products ED multiplied by the SCET quark spinors and gluon polarization 
vector. As already mentioned the tree matrix element of the “right insertion” of Qi 
vanishes due to colour, so A- 0) = 0 for Q\. In this case the 1-loop amplitudes are 
IR-finite and can be evaluated in d — 4 (after UV renormalization is applied). Hence the 
evanescent terms i — 2,3,4 vanish. The “right insertion” of Q 2 has A[ 0 ^ ^ 0 for i — 1,2 
[see mi and the 1-loop amplitudes are IR-divergent. The i = 1,2 1-loop terms have a 
1 /e 2 singularity, proportional to the tree matrix element, as follows from the universality 
of soft singularities. A^ 1 ' 1 has a 1/e pole, while A{ 1) turns out to be IR-finite. The IR 
divergences cancel when the QCD amplitude is related to the matching coefficient H u 
through (ED as explained in the following. 


3.3 IR subtractions 

We start from the matching equation (0) extended to include the evanescent operators 

Q = T l *O l + Y. H ?*°f- ( 22 ) 

i 

Convolutions, which may involve one or two integrations, are now represented by an 
asterisk. Since we work with matrix elements in states with definite momentum it is 
convenient to use the momentum-space representation. Expanding all quantities to the 
1 -loop order, making use of m and Of = Oj, we obtain 

Y, (4 1} + A'f) <Oi) (0) = T l U(O l )W + T m * (OT ] 

i 

+ Y + Hj m * (OYA (23) 
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The factorizable contribution on the left-hand side comes from 1-loop diagrams with 
no gluon lines connecting the XX part of the diagram to the part. It is canceled 

by the term T I(0 ) * (O 1 )*- 1 ) on the right-hand side, since the 1 -loop matrix element of 
0\ contains exactly these diagrams in the coefficient function C^ 1 ' 1 in its definition (J3J). 
The UV-renormalized 1-loop matrix elements of the O, are given by 

<0,) (1 > = V (M« R + Z\f) (24) 

3 

/i\n 

where is the bare matrix element, which depends on the IR regularization scheme 

R , and zjP the matrix kernel of ultraviolet renormalization factors. When dimensional 
regularization is used for UV and IR singularities as was done in the calculation of A^\ 
the bare matrix elements vanish, since the 1-loop diagrams are scaleless. Hence, inserting 
m into (1231) . using (J3J) and Hj 1 ^ = Af \ we obtain 

4 1} = H{ m + J2 A i 0) * z a + (-2)r I(1) C'^ 1)(0) (25) 

i=1 

by comparing the coefficient of (Oi)®. We also used that A\ 0 ^ is zero for i — 3,4. The 
renormalization constants for the evanescent operators are determined by requiring that 
the IR-Hnite matrix elements (Oj) (i = 2,3,4) vanish jT3J [o'. Here “IR-Hnite” means 
the matrix element computed with any IR regularization R = off other than dimensional 
and with dimensional regularization applied only to the UV singularities. According to 
(12411 this fixes = — M$ oS . Hence the 1-loop short-distance coefficient of the physical 
(non-evanescent) operator O il is given by 

H n{i) = A A) _ A (°) * + 4°) * M W° S + 2 T i(i) c (Bim _ ^ 

Note that since the IR-hnite matrix elements of the evanescent operators have been 
made to vanish, only the term i = 1 survives in (I 22 [) . It is therefore not necessary to 
determine the coefficient functions H\ l for i ^ 1. We also note that the renormalization 
constant Z^\ is hnite, and that is independent of the apparently arbitrary IR 

regulator. This is because the mixing of an evanescent operator into a physical operator 
arises through the multiplication of an ultraviolet 1 /e pole with a term of order e from 
the Dirac algebra, both of which are independent of the IR regularization. The 1/e 2 
poles do not contribute to operator mixing due to their universality. 

Eq. (12611 provides the final result for the two of the four matching coefficients asso¬ 
ciated with the right insertions. We briefly discuss the subtraction terms in (HU). First 
note that for the right insertion of Qi the tree amplitudes A® vanish, hence (Infill is 
simply H ll[ ' l) = A^. This is consistent, since for this case A^ is IR-hnite as observed 
above. For the right insertion of Q 2 , all three subtraction terms in (HU) are present. 
Due to factorization in SCET, the renormalization of an operator [xx] [£ A ±cih v \ falls 
apart into a renormalization factor for the collinear-2 bracket [xx] and one for a B 1-type 
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Figure 4: SCET 1-loop diagram contributing to M 2 i^° ff , the mixing of 0 2 into 0\. 


current. z[\ } is therefore determined by the requirement that the light-cone distribution 
amplitude of M 2 and the jet function are defined in the MS scheme. This gives Z as 
the product of the Brodsky-Lepage kernel m and the renormalization kernel Z\\ for the 
Bl-type current (first paper of 0, [ID])- Subtracting from removes the 

IR singularities such that H II( k is finite as must be for a short-distance coefficient. The 
third term on the right-hand side involves the computation of the 1-loop matrix element 
of the evanescent operator 0 2 . We fold that only a single diagram, shown in Figure 01 
contributes to M 2 ^ off , such that M^ oS is proportional to the spin-dependent part of the 
Brodsky-Lepage kernel. Explicitly, 



M. 


(l)off 


21 


a s Cp 

47T 



(- 1 ) 

N c u'u' 


(- 8 ) 



u 


u’) + - 6(u’ 
u 



a s Cp 

A^W C 


(- 8 ) 



In u 
u 


(27) 


The fourth term 2 follows from = —1 lT) and [I] 

T I(1) = y ^ (28) 

4vr 1V C v ’ 

with V(u) given in (PHI) below. Note that H u ^\u,v) is a function of two variables, but 
like the tree contribution the two subtraction terms (EH), EED depend on u, but not on 
the momentum fraction v related to the collinear-1 momenta. 


3.4 “Wrong insertion” 

The other two matching coefficients are related to the wrong insertions of Q 12 as in the 
right diagram of Figure 0 We would like to express them as the coefficients of the same 
SCETj operator O u (t, s), but the QCD calculation involves Dirac matrix products with 
a different contraction of spinor indices corresponding to 


O u {t,s) 


— [(£W cl )(tn_) W cl ] (sn + )(1 - 7s ) (W c k) 
7 ±(! — 75) ® (I + 75) 


( X W c2 )(1 + 75 )^ 


(29) 
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In the second line we introduced again a short-hand notation that highlights the Dirac 
structure. The symbol 0 means that the spinor indices are contracted as in [£y] [xh v \. 
We deal with the required Fierz transformation and evanescent operators simultaneously 
by introducing the operators 

Oo = ^(1 - 75 ) ® ^y(l - 7 s)7jS 
di = t±(i — 75) ® (1 + 75), 

C>2 = 7 ±(! - 75) ® (1 + 75 ) 7 ± 7 ±a, 

O 3 = 7jl7±7±(! - 75) <8» (1 + 75)7±«7±/3, 

O 4 = 7±7±7±(1 - 75) ® (1 + 7s)7±7±«7±/37±7- ( 30 ) 

Here O 0 is the short-hand for 0 11 . The basis is chosen such that Oi and Oo are Fierz- 
equivalent and O 2-4 vanish in four dimensions. Hence we have one physical operator, 
O 0 , and four evanescent operators, 0 \ — O 0 , and 02 - 4 - The tree matrix element is now 
given by 

«l)n, = T (1(0,) + i{0 2 )) , (31) 

N c \u u ) 

and the 1 -loop amplitude reads 

{<?>„, = £ (if + if) (O^K (32) 

1=1 

This does not contain (O 0 )(°), since all diagrams have the “wrong” Fierz-ordering. Pro¬ 
ceeding as before and requiring that the infrared-finite matrix elements of the four evanes¬ 
cent operators vanish, we find that (OKI) is replaced by 

H m) = £0) _ i[°) * + 4 0) * m£ ) oS + iS 0) * (. m[\ ) oS - MoQ )off ) + 2 T^ l) Cf + im (33) 

with M 2 ( ! )off the bare 1-loop mixing of O 2 into Oi, and = Z,The new term 
*( k M[^ ) oS — Mqq° S ) involves the difference of the mixing of 0\ and Oo into themselves. 
This difference is finite and independent of the IR regulator for the same reason that 
is. There is one subtle aspect hidden in © that requires explanation. As in 
© we would like to cancel the factorizable QCD diagrams against the matrix element 
of O 1 , but the two terms appear in different Fierz-orderings. The consequence of this is 
that there should be an extra term related to the factorizable diagrams on the right-hand 
side of ©. Using that at tree-level only , A^f ], A^j are non-zero, it is given by 



However, we find that this term vanishes, hence © is correct. The subtractions — Aj° * 
and 2 T i ( 1 )C'| B1 ^°) are identical to the corresponding terms for the right insertion. 
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The other two terms are once again related to an integral over the spin-dependent part 
of the Brodsky-Lepage kernel. Explicitly, they read 



(35) 


Despite the different Dirac algebra and subtraction structure we find that the final result 


for the matching coefficient related to the wrong insertion of Qi (Q 2 ) is identical to the 


one for the right insertion of Q 2 (Q i). 

3.5 Matching coefficients 

Here we give the final results for the matching coefficients (hard spectator-scattering 
kernels) H n (u,v) in the convolutions (fTj). (fTTHh The kernels are independent of the 
mesons but dependent on their flavour quantum numbers. To make them explicit, we 
write {[q s QM 1 ][QM 2 < lM 2 ]\Q\[Qsb]) for {MiM 2 \Q\B) to indicate flavour. Then 

- for ([q s D][uu}\Qi\[q s b]), the contribution of Q i to the colour-suppressed tree ampli¬ 
tude a 2 (MiM 2 ) , and for ([q s u][uD]\Q 2 \[q s b]), the contribution of Q 2 to the colour- 
allowed tree amplitude cki(M 1 M 2 ), we have 



(36) 


- for {[q s u][uD]\Qi\[q s b]), the contribution of Q i to the colour-allowed tree ampli¬ 
tude cki(MiM 2 ), and for ([q s D][uu]\Q 2 \[q s b]), the contribution of Q 2 to the colour- 
suppressed tree amplitude a 2 (MiM 2 ): 



(37) 


Here ri(u, v), r 2 (u,v) are given by 
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where we defined 


F(v, w) = 2 Li 2 ( -— + 2 Li 2 (u>) — Li 2 (™) H— In 2 — + in In —, 

V w J 2 v v 

V(u) = 6 In — 18 + 3 In u — in\ 

H z \ u J 

( 2 In 'll 1 

+ |2 Li 2 (w) — In 2 u A -—-(3 + 2m) In it — (u u)>. (40) 

The expressions for rq and r 2 constitute the main technical results of this paper. In 
applications the kernels always appear in convolutions. In the following we perform the 
convolution integrals analytically and obtain a compact representation for the (topolog¬ 
ical) tree amplitudes. 


4 Tree amplitudes with NLO spectator scattering 


The complete 1-loop correction to spectator scattering is the convolution of the hard- 
scattering kernels H\ ] 2 with the jet function 

47t a s Cp 1 


Jii ( z ; v , u ) = 


N c rriBUJV 


CL 


S(z-v) + 

47T 


(41) 


Inserting this and (ESI), EH) into m , and expanding to order a 2 , we obtain 


Ti (a ;,u,v) 


'KOLgCp 1 

N% rribUJV 


1 f 1 

ri(u,v) + - dz j\\ (z; v, u) 
u Jo 

r 2 (u,v) (i = 2) 


Ots_ 
U 47T 
OL s 

An 


(■ = i) 

(42) 


(The integral /q dzj\\ (z\ v, a;) is given analytically in appendix B.l of PHI-) The spectator¬ 
scattering contribution to the tree decay amplitudes ai i 2 (MiM 2 ) in the standard nor¬ 
malization is 


Sl = 7^%) l Jo dudvT i l ^i U i V )^B+(u)(t)MA v )^M 2 (u), (43) 

see ©- More precisely, accounting for the Wilson coefficients in the effective weak 
Hamiltonian El, we have ai(M 1 M 2 )| S p — Cj S 2 + C 2 Si, a 2 (M 1 M 2 )| sp — CjSj + C 2 S 2 . 


4.1 Expansion of convolutions in Gegenbauer moments 


The integration over the spectator quark momentum fraction u is simple, because u 
appears only in the jet function, or as the over-all factor 1/a;. The dependence on the 
light-cone distribution amplitude of the B meson is encoded in the inverse moment 


1 



Jo u; 


(44) 
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and the logarithmic moments 


(L"> 



(45) 


up to n — 2. The light-cone distribution amplitude of a light meson, <p M , is conventionally 
expanded into the eigenfunctions of the 1-loop renormalization kernel, 




6xx 


l + Y.<CW\2x-l) 

n= 1 


(46) 


where a™ and are the Gegenbauer moments and polynomials, respectively. The 

integrals over u and v can then be performed, the result being represented as a double 
expansion in the Gegenbauer moments of Mi and M 2 . Often, there appear the quantities 


Am = 



4>m(x) 

3x 


1 + £(-!)” “f. 

n= 1 


A M = 



3x 


i + E 


M 


n=1 


(47) 


We give the final result for the tree amplitude parameters aq i2 (MiM 2 ) in the notation 
of [3] [eq. (35)], including the 1-loop vertex correction 14 (M 2 ) not related to spectator 
scattering, 


OiiMiMz) = Ci 


t ^ -I ■ i C,+i gl q (4 f . . 


N r 


N r 4vr 


TlOigC f 

~nT 


9/mi/b 


Q±i hi(M l M 2 ) + Ci h 2 (M 1 M 2 ) 


(48) 


The upper signs apply when i is odd (here simply i = 1), the lower ones when i is even 
(here i = 2). The spectator-scattering mechanism is encoded in the two objects 


hi{MiM 2 ) 


h 2 {M\M 2 ) 



j R 1 (M 1 M 2 ) + Am 2 J(M 1 ) , 


cv 


(49) 


such that the a s terms in these expressions extend the result given in PI, and r™ 1 A Ms X H 
denotes a power correction included in the definition of Hi(M\M 2 ) in Pj. Performing the 
convolution integrals in a double Gegenbauer expansion as described above, the hard- 
collinear 1-loop (o^) correction is given up to the second Gegenbauer moment in terms 
of (|H3, appendix B.l) 


J(Mi) = ^ f — <j> Ml (v) j — 0b+( u) j drj\\ (t;u, 
3 Jo V Jo id Jo 


id 


= V) + (A? + q (L) + "»- ^-4(3) 

3 X 1 \ 3 9 / W 18 9 3 V ; 
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(50) 
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(Here we have set nj = 4, Tj = 1/2 and Cj? = 4/3 and Ca = 3.) The new hard correction 

is in 

1 r 1 . . , , , , , , r k (u,v) 


Rk(M 1 M 2 ) = - dudv (j) Ml (v) 4 >m 2 ( u ) 
9 Jo 


(51) 


Integrating the kernels (13)811 . (I, ‘-SOI) , and truncating the Gegenbauer expansions after n — 2, 
we obtain 
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Here we defined i = In(m^//r 2 ). The finiteness of R\^ proves factorization of spectator 
scattering to the 1-loop order. For /i = rrib (£ = 0) the magnitude of the correction 
and importance of the higher Gegenbauer moments can be seen from the numerical 
expressions 

R\ = 0.6047 + 10.9683 i + (25.36 + 23.19 i) of 2 + (46.39 + 29.41 i) af 2 
+ (-6.10 + 14.08 i) af 1 + (-12.84 + 15.08 i) af 1 
+ (22.69 +23.81 i)af 1 af 2 + (44.63 + 30.83 i) af^f 2 
+ (19.21 + 23.80 i)af 1 af 2 + (43.62 + 31.35 i) a^af 2 , (54) 

R 2 = -8.2259 - 5.2906 i — (4.31 + 8.43 i) af 2 - (2.59 + 8.24 i) af 2 
-(7.16 + 6.58 i) af 1 - (7.83 + 7.04 i) af 1 
— (2.36 + 10.95 i) af x af 2 - (0.83 + 11.28i) 2 
-(3.63 + 11.98 i)af x af 2 - (2.37 + 12.60i) a^af 2 . (55) 

For phenomenology the most important 1-loop spectator-scattering correction is the 
term C\R\ involving the large Wilson coefficient C\ in the expression for a 2 - From 
the above expression for R\ we learn that this contribution has a large imaginary part, 
while the real part seems to be accidentally small. The higher Gegenbauer moments 
have relatively large coefficients. Roughly, the magnitude of the correction to h 1 (M 1 Af 2 ) 
is 14a s /(47r) « a s ~ 0.3, i.e. a 30% correction relative to the tree approximation. A 
detailed numerical analysis of the tree amplitudes will be performed below. 

4.2 Scale issues 

Up to this point we did not make explicit the scale dependences of coupling constants and 
parameters. The tree amplitudes ol\^ 2 themselves are scale-independent, but the Wilson 
coefficients Ci, strong coupling a s , static B decay constant /#, light-cone distribution 
amplitudes of the light mesons (hence the Gegenbauer moments af 4 ), as well as the B 
meson distribution amplitude moments A b, ( L n ) are all scale-dependent. Expressing all 
these quantities at the scale /i equal to the one that appears explicitly in the 1 -loop result 
for the hard-scattering kernels is a legitimate choice. However, with any single scale one 
or another kernel will contain parametrically large logarithms lnmfc/AQCD- 

For the following discussion we assume that the hadronic quantities af 1 , Ab as well 
as the logarithmic moments 

<+(/+) = A B (/i 0 ) / — (p B +(u) ln n — (56) 

JO W LU 

related to ( L n ) are given at a reference scale /i o of order of the hard-collinear scale 
(m{,AQCD) 1//2 - The first line of f|48|h which corresponds to the form factor term in the 
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factorization formula CD or (JHD is easy, since it contains only a hard correction. There are 
no large logarithms when the Wilson coefficients and vertex kernels Vi(M 2 ) are evaluated 
at a common scale / 1 & of order m b . Since the tree approximation is independent of the 
Gegenbauer moments, the Gegenbauer moments in Vi(M 2 ) are evolved from /j, 0 to /ib 
in the leading-logarithmic (LL) approximation with the 1-loop anomalous dimension 
matrix. 

The spectator scattering term is more involved. In order to resum large logarithms 
one should perform the substitution 

GO) H n (fj.) * * [f B (j) B+ }(fi) * <j> Ml (fj) * </>m 2 O) 

=>• C(Hb) H U (fi b ) * U{lXb,fl he) * J\\ (h-hc) * [/n0B+](hhe) * 0Mi(hhc) * 4>M 2 {d he) (57) 

where /Zh c ~ (mfoAQCD) 1 / 2 is of order of the hard-collinear scale, and f/(/z&,/Zh c ) is the 
evolution function for the SCETi operator O n . While in the first line either H n (/j.) or 
Jm(/z) contains large logarithms, neither nor J\\(/J<hc) does. The evolution function 

for O 11 factorizes into Ubl related to the Brodsky-Lepage kernel in the collinear-2 sector 
and the evolution function Un for a Bl-type current. Since Ubl(/W /-the) * </>M 2 (hhe) = 
<Am 2 (/•*!>)) we can rewrite the previous expression as 

C(fJ. b ) * U\\(nb, /Zhc) * J||(//hc) * [Mb +1 (Mhc) * 0 Mi(/rhc) * 0 M 2 (hfo)- ( 58 ) 


In this expression the Gegenbauer moments of M 2 must be evolved from /z 0 to fib in the 
next-to-leading-logarithmic (NLL) approximation using the 2-loop anomalous dimension 
matrix, since the Gegenbauer moments appear already in the tree approximation to 
spectator scattering. We have implemented the 2-loop evolution based on the results 
of JT7J. The evolution preserves the truncation of the Gegenbauer expansion due to the 
triangular structure of the anomalous dimension matrix. Up to the second Gegenbauer 
moment the evolution is obtained from 


d 

d ln/z 
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with (putting rif — 4) 
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(60) 


The Gegenbauer moments and B meson distribution amplitude moment A b in J ||(/Zh c ) * 
[/u0s+](hhe) * 0Mi(hhc) can be obtained from the input values at /z 0 by a fixed-order 
1 -loop relation, because no large logarithms appear when these factors are evaluated at 
/z hc . Ib is obtained from the physical decay constant fs by a HQET conversion factor 
using the 2-loop approximation for the anomalous dimension and the 1-loop matching 
coefficient, since the matching to fs is done at the large scale /z&. To complete the 
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Parameter 

Value/Range 

Parameter 

Value/Range 

A MS(5) 

0.225 

l^b 

°0 ^ 
^ CN 
+ 1 
OO 

m c 

1.3 ±0.2 


1.5 ±0.6 

m s (2 GeV) 

0.09 ±0.02 

fs d 

0.20 ±0.03 

m b 

4.8 

/+*( o) 

0.28 ±0.05 

rh b (rh b ) 

4.2 

A s (l GeV) 

0.35 ±0.15 

\V A \ 

0.0415 ±0.0010 

ui(l GeV) 

1.5 ± 1 

\v ub /v cb \ 

0.09 ±0.02 

<r 2 ( 1 GeV) 

3 ± 2 

7 

(70 ± 20)° 

a£( 2 GeV) 

0.1 ±0.2 


Table 1: List of input parameters. Dimensionful parameters are given in units of 1 GeV. 


evaluation of (1551) one would now require the 2-loop anomalous dimension kernel for the 
Bl-type current to evaluate U\\(fi b , /-the) hi the next-to-leading-logarithmic approximation. 
Only the 1-loop anomalous dimension is available (first paper of [TQ]). We therefore 
implement an approximate procedure analogous to that in EH for heavy-to-light form 
factors. We evaluate i7 n (/i) at the hard-collinear scale /ihc except for the terms involving 
logarithms related to the scale dependence of the Wilson coefficient C and of 0 m 2 > which 
remain at /i b . To this expression we add the series of all leading logarithms summed to 
all orders in perturbation theory omitting the terms already included in L/ n (/i hc ). The 
structure of these terms is identical to eq. (117) of H3J- 

Finally, we match a s from a 5-flavour to a 4-flavour theory at the scale [i b . Quantities 
evaluated at /i hc are computed in the 4-flavour theory, quantities at /if, in the 5-flavour 
theory. 

5 Numerical analysis of a\ and oq and application to 

the B —> 7777 system 

5.1 Input parameters 

For our numerical study of the tree amplitudes cq and «2 we employ the input parameters 
shown in Table[T| With respect to (3] we have updated the value of |W&|, and reduced the 
error estimate of aj. The heavy-quark masses m b , m c are now interpreted as pole masses 
with m b = 4.8 GeV. The list is extended by the logarithmic moments of the B meson 
distribution amplitudes, for which we use values in the ranges obtained from QCD sum 
rules or models for the shape of the distribution amplitude |^|. The hard-collinear scale 
/ihc and the hard matching scale /if, are varied independently within the given ranges. 
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5.2 Tree amplitudes aq and 

Numerically, we obtain the tree amplitudes 3 * * * 

ai(7T7r) = 1.015+[0.025 +0.012* ] v - (^ 2 ) {[0-014 + 0.0245 a2 ] LO 

+ [0.024 + 0.0445 a2 + 0.0202 + 0.050<5 a2 ] NLO + [0.009 + 0.007<5 a2 ] tw3 } 

= 0.99218;™ +(-0.00718;“®) i, (61) 

a 2 (7T7r) = 0.184 - [0.152 + 0.0772 ] v + {[0-088 + 0.1505 a2 ] LO 

+ [0.029 + 0.130h a2 + 0.0342 + 0.100<5 a2 ] NLO + [0.056 + 0.04l5 a2 ] tw3 } 

= 0.20518:^ +(-0.04318:81!) 2'- (62) 

In these expressions we separated the tree (first number), vertex correction (indexed by 
V) and the spectator-scattering correction (remainder). The latter is further divided 
into the tree (ay, indexed LO), 1-loop (a^, indexed NLO), and twist-3 power correction. 
The 1-loop correction is the sum of the jet function and hard correction, see (USD- We 
also made explicit the parameter dependences that are responsible for the bulk of the 
theoretical uncertainty given in the last line in the expressions for ay 2 . (Theoretical er¬ 
rors computed from the ranges in Tableware added in quadrature.) The most important 
such parameter is the combination 


_ 9/tt/b 

Vsp ~ m b f^(0)X B 7 

which normalizes the spectator-scattering term as can be seen from (1481) . The second 
most important parameter is the second Gegenbauer moment of the pion distribution 
amplitude. This dependence is shown (for spectator scattering only, where it is impor¬ 
tant) to linear order in the deviation 5 a2 = aj(2 GeV) — 0.1 from the default value. The 
result for the two tree amplitudes is shown in Figure O which also displays a comparison 
with the leading-order and next-to-leading order ( 1 -loop vertex correction, tree spectator 
scattering) result, as well as the main parameter dependences. It is evident from (liHT) . 
© or the Figure that the 1-loop spectator-scattering correction is significant, but not 
large enough to put the perturbative approach into question. For a 3 the 1-loop correc¬ 
tion exceeds the tree spectator correction, because the 1 -loop correction is multiplied by 
the large Wilson coefficient Ci. However, the correction is small in absolute terms. For 
a 2 the correction amounts to about 35% of the tree spectator-scattering term. Since the 

3 The following numbers differ from m , since in m the hard 1-loop spectator-scattering correction 

has been added to the program used in Pj to allow for a direct comparison with the scenarios defined 

there. For the present numerical evaluation the program has been substantially changed in order to im¬ 

plement the various scale dependences as described in Section 1Q1 In particular, the Wilson coefficients 

Ci are now evaluated at /ib, and always in the NLL approximation. 
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Figure 5: The tree amplitudes aq(7r7r) and a 2 ( 7r7r ) represented in the complex plane. 
The dark (black) diamonds show the LO, NLO, and partial NNLO approximations. 
The latter includes the new 1-loop correction to spectator scattering and is shown with 
error bars. The dark square represents the parameter set ‘G’, which provides a good 
description of the experimental data on branching fractions as discussed in Section 15.hi 
The grey (blue) triangles show the variation of the tree amplitudes, when A b takes the 
values 0.2 GeV to 0.5 GeV in steps of 75 MeV, such that the triangles in the direction of 
the point ‘G’ correspond to smaller values of A#. From each triangle emanates a set of 
grey (red) points that correspond to varying a 2 from —0.1 to 0.3 in steps of 0.1 for the 
given value of A b- Here points lying towards ‘G’ correspond to larger aj. 

imaginary part is generated only at NLO, it is best compared to the imaginary part of 
the vertex correction V. This shows that the spectator-scattering correction at order 
is almost as large as the vertex correction at order a s , but comes with an opposite sign 
such that the phases tend to cancel. 

With the perturbative approach thus validated through the size of the 1-loop correc¬ 
tion, it is evident from the Figure that the dominant uncertainties are due to hadronic 
input parameters. The uncertainties in /#, A b and f+ n (0) do not exclude that r sp is 
a factor of 2 larger than its default value 0.412. In fact, it appears that the data on 
B —> 7T7T branching fractions require such an enhancement [3|. Until some of these pa¬ 
rameters are better determined (from theory, from other data, from fits to non-leptonic 
data) there remains a large uncertainty in the colour-suppressed tree amplitude a 2 . The 
colour-allowed tree amplitude, however, is predicted to be close to 1 with an uncertainty 
of 10% even with present parameter inaccuracies. 
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5.3 B —» 7T7T branching fractions 

We confront our new (partial) NNLO results with the experimental data on the three 
tree-dominated B —> ixix branching ratios. The B —> nn amplitudes are given by 


V ^2 j\q —>7r 71*0 
—>7r"*" 7T 

A§°^7r°7r 0 


i + 02 

i % m l/»/b(°) {uiKdbi + a; 

i 0= m !/*/f'( o ) {ViiKi[aj - a? 


+ U,V>S}, 
- KiK>s}, 


(64) 


not showing some smaller amplitudes that are taken into account in the numerical evalu¬ 
ation of the branching fractions below. The theoretical computation includes the 1-loop 
correction to spectator scattering in the tree amplitudes, 0 : 1 , 2 , but not to the QCD 
penguin amplitudes a 4 ,c . For these (and the smaller amplitudes not shown) the NLO 
expression ma updated to include the scale-dependent parameters aj and fs in the 
LO approximation for spectator scattering is used. 

The standard input parameter set does not provide an adequate description of B —> 
7T7T data. Rather the data favours a smaller value of |14&| which reduces the overall 

normalization of the amplitudes, and a significantly larger contribution from spectator 
scattering, which increases oq (see Figure 0) (3J • We find that the parameter choice ‘G’ 
with 


|K&| f+B 0) = 8.10 • 10- 4 = 0.775 \\ v ub \ 0) 


. def ’ 


r sp = 1.96 \r 


spjdef 


(65) 


and 0^(2 GeV) = 0.3 yields a good description of data. The required parameter modifi¬ 
cation is most likely related to a smaller value for the B —> n form factor and a smaller 
value of As, but other small modifications may add up to the combined effect. The 
new parameter selection ‘G’ is similar to scenario S4 defined in [3J, and falls within the 
ranges for the individual parameters specified in Table El With (1H5I) we calculate the 
CP-averaged branching fractions 


10 6 Br(R" -> 7r _ 7r°) 
10 6 Br(R° -»■ vr+Tr-) 
10 6 Br(R° -> 7r°7r°) 


^■5io!3(CKM)12!4(hadr.)12;I(pow-) [5.5 ± 0.6], 

5.0+°;®(CKM)+^(hadr.)+^(pow.) [5.0 ±0.4], (66) 

0.73i^I(CKM)l°;“(hadr.)^;l(pow.) [1.45 ± 0.29] 


with the experimental averages reproduced in brackets EDI The corresponding tree 
amplitudes aq = 0.92 — 0.05i and ct 2 = 0.51 ± 0.03i are shown by the points marked 
‘G’ in Figure El which implies that the ratio of the colour-suppressed to colour-allowed 
amplitude C/T = a^/aq = 0.55 ± 0.07i is large. By construction the branching frac¬ 
tions with charged pions in the final state are in excellent agreement with data. The 
B —> 7r°7r° branching fraction is still somewhat low, but the theoretical uncertainty is 
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large. In computing the theoretical errors we did not include here the uncertainties in 
\Vub\, /+ 7 r ( 0 ), /b, A b and aj, because these 5 input parameters appear only through 
m- The remaining uncertainties are divided into groups from \V c b\, 7 (CKM); the 
renormalization scales fib, and quark masses m c , m c (hadr.); and X H , Xa (weak 
annihilation) parameterizing non-factorizable power corrections (pow.). The decays to 
the final states 7 r + 7 r“ and 7 r° 7 r° are sensitive to 7 as can be seen from the CKM uncer¬ 
tainty. The dominant errors come from the hard-collinear factorization scale fi^c, and 
from power corrections. We postpone a detailed assessment of the theoretical status after 
the calculation of the 1 -loop spectator-scattering correction to the penguin amplitudes, 
which may be important for B° —» 7 r° 7 r°. We also expect the spectator-scattering phase 
in the penguin amplitudes to affect the direct CP asymmetries, and therefore do not 
discuss them now . 4 


6 Conclusion 

We computed the 1-loop hard spectator-scattering correction to the topological tree 
amplitudes in non-leptouic B decays by matching the current-current operators Qi ,2 
from the effective weak interaction Hamiltonian to the relevant operators in SCETi. To¬ 
gether with the 1 -loop calculation of the hard-collinear jet function [3 HU], the spectator¬ 
scattering contribution is now complete at order a also providing the first and perhaps 
dominant contribution to a full NNLO computation of the decay amplitudes. Unless the 
1-loop term is enhanced by a large Wilson coefficient, it is of order (30 — 50)% relative 
to the tree term, depending on moments of light-cone distribution amplitudes. This 
yields a visible enhancement of the spectator-scattering amplitude without changing the 
qualitative features of the previous approximation. 

The very fact that the perturbative correction can be computed, and that the expan¬ 
sion seems to be reasonably behaved, is significant, since it shows that a) factorization 
holds theoretically, i.e. IR singularities factorize as predicted and convolution integrals 
converge, and that b) perturbative expansions of the spectator-scattering contribution 
are under control, an issue that has at times been a point of controversy (second paper of 
mi. eat Our calculation therefore shows that the theoretical accuracy is now limited 
by the uncertainties of input parameters for the factorization formula, and encourages 
efforts to determine better key hadronic quantities such as form factors, and moments 
of light-cone distribution amplitudes. 
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